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IL. INTRODUCTION 


Atmospheric turbulence produces optical phase and amplitude fluctuations that se- 
riouslv degrade laser beam propagation and the performance of imaging svstems. In 
order to perform svstems analvses of optical propagation, the spatial and temporal dis- 
tribution of atmospheric optical turbulence must be understood. Over the past two and 
a half vears optical turbulence measurements in the form of isoplanatic angles, 0,, and 
transverse coherence lengths, rj. have been made at a variety of sites in the continental 
United States, as well as Maur, Пахжап апд Vieques, Puerto Rico. Though a dominant 
application of these measurements involve laser technology and astronomy, the inter- 
pretation of the results are very valuable to the meteorologist. The optical data repres- 
ents an instantaneous summary of the optical turbulence along a vertical path. 
Imbedded within this vertical cross-section is information describing various synoptic 
and mesoscale phenomena such as the jet stream, gravity waves, low level jets and 
boundary laver circulation. By correlating meteorological measurements and optical 
turbulence parameters, the practical cause and effect relationships can be assrmilated 
into prognostic tools. 

lx nied States Air Force Relay Mirror Experiment (RME) requires a critical 
evaluation of atmospheric turbulence over the Maur. Hawail experiment site. RME is 
designed to "assess how accurately an orbiting mirror can place (a) reflected (laser) beam 
EE Uround target with low beam wander" (Nelson, 1988). Observed beam jitter and 
Wander are a product of both the laser svstem and the atmospheric turbulence. Accurate 
measurements of the atmospheric contributions are vital for evaluating the optical sys- 
ливре огтапсе. 

Between August 1987 and July 1988, the Naval Postgraduate School (NPS) per- 
formed six measurement sessions of approximately one week duration each in order to 
quantify the RME target board site’s optical turbulence characteristics. NPS optical 
instruments were used to acquire atmospheric coherence length, 7), and 1soplanatic angle, 
0,. measurements near sea level at Kihei, laut. In addition to meteorological data from 
the National Weather Service Rawinsonde Observation Stations. synoptic charts and the 
GOES-WEST Infrared Satellite Images, four of the six Kihet data sets had on-site 


rawmsonde measurements. 


This thesis attempts to correlate the RME meteorological measurements with the 
optical data. Initially, the National Weather Service measurements, charts and satellite 
images Were used to compare meso- and synoptic scale events with the optical results. 
This method provided qualitative results. 

The major thrust of this thesis, however, materialized with the availability of high 
resolution rawinsonde soundings. These profiles finally provided the resolution neces- 
sary for attempting the direct computation of optical parameters from the rawinsonde 
data. The August and December 1987 sessions used the VIZ WO-8000RP + rawinsonde 
svstem. With a vertical resolution of approximately 35 m, measurements were recorded 
upto 12 km. The April and July 1988 sessions made use of the Vaisala DigiCORA MW 
l] system. This system profiled meteorological data up to approximately 24 km with а 
sampling rate of 2 seconds (surface to 2 km) and 10 seconds (above 2 km). The $ m and 
40 m vertical resolutions of the DigiCORA system Were considered reasonable for com- 
puting the optical turbulence directly from the rawinsonde data. 

The optical variables measured. r, and 0,, depend on an integral of the index of re- 
fraction structure parameter. C2, along the vertical path. The parameter Cass s 
tional to the temperature structure parameter C} which can be computed from the 
rawinsonde data. Four methods for computing the temperature structure parameter C} 
from the rawinsonde data were investigated: 


1. the Simple Potential Temperature Differences (SPTD) Equation, 


tJ 


the Potential Temperature Diflerences with a Running Mean. 


"12 


the Bulk Gradient Equation (BG) and 


4. conditional data sampling based on the Richardson Number. 


Lach Ci method was tested by calculating л, and 0, from the C; profilce оси ин 


were ultimately compared to the optical measurements simultaneousiv acquired over the 
} ] ] 


Same site. 


tl 





H. BACKGROUND 


A. HISTORICAL REVIEW OF CORRELATION RESEARCH 

International scientific interest in the meteorological influence on "seeing, a rating 
system based on the angular size of a star image, flourished in the mid to late 50's. 
Following a brief decline, interest was reestablished in 1962. The new upper atmospheric 
and outer space environmental data introduced by artificial satellites and space probes 
provided a catalyst for the development of atmospheric physics. One of the dominant 
questions addressed by the melange of scientific backgrounds required in the atmo- 
spheric physics field involved how meteorological conditions affect “seeing” (Kucherov, 
1966a). 

In October 1962, the issue of observatory site selection and related astroclimatic 
topics were discussed at the International Astrononucal Union Svmposium held in 
Rome. Geophvsictsts, meteorologists and astronomers worked in cooperation to study 
the interrelationships between atmospheric conditions and astrononucal “seeing”. Par- 
ticipants resolved that placing telescopes at certain heights above the ground would 1m- 
prove image quality (Kucherov, 1966a). The following summarizes some of the 
pertinent articles presented at the Third All-Union Conference held four months earlier 
in Kiev. 

Kucherov (19665) observed "shimmer angle" and meteorological conditions at one 
seashore and two inland sites over a period of approximately two уеагѕ (1960-1962). 
The results of that study confirmed the relationship between synoptic fronts (mostly cold 
fronts) and astronomical “seeing”. It should be noted that 25% of the cases display no 
relationship. 

Bol’shakova and Darchiva (1966) investigated the causes of “seeing” condition 
changes. While thev acknowledge the synoptic fronts as an influence on “Seeing”, they 
suggest that “seeing” is affected by frontal passages only when air masses with optical 
characteristics differing from the original advect over the site. Using statistics, thev 
support their claim. To a first approximation, 5075 of the cases show the synoptic front 
causing a deterioration in “seeing”, or alternatively, absence of a synoptic front ensures 
high “seeing” conditions. The rest of the cases display the reverse situation. In their 


opinion. the analvsis must be repeated taking into account the type of front, the incli- 


nation of the frontal plane, the direction of the front’s motion, and ideally, other detailed 
meteorological information defined in various atmospheric lavers. 


Zinchenko (1966) comes to the following conclusions: 


1. The lower-atmospheric turbulence is of decisive significance in all that concerns 
“secius 


a. "optical inhomogeneities are most effective in the ground laver, where the air 
density is highest:” 


b. "turbulence intensitv is greatest in this layer апа 
c. “eddies of minimum size occur in the ground laver". 

2. Contradictory results exist with the applicability of Richardson number(Ri) as a 
criterion of high atmospheric turbulence. 

The main reason supporting the latter point 15 that turbulence intensity is not a 
single-valued function of Ri, though a dependency on Ri exists. Zinchenko explains, 
When Ri< 1, one can qualitatively accept that turbulence set in: and, the amount of ac- 
cumulated eddy enerev in steady-state conditions increases with the decrease іп Кі. 
(Note: Current understanding requires Ri< 0.25 for the above to be true.) опен ин 
oi data prohibited an exact definition of the laver directly responsible for the blurring 
of the stellar diffracuon disks. Zinchenko found that turbulence does show a tendency 


to be restricted to the 0 5-2. 0 bn F o. 


D. THE RELAY NIIRROR ENPERINIENT 

The Relay Mirror Experiment (RME) is a United States Air Forces s мн ни 
monitored by the Kirtland Air Force Weapons Laboratory. Albuquerque, New Mexico. 
RME is designed to “assess the ability to place a laser beam on a ground target with low 
beam wander” (Nelson. 1988). As described by Nelson. RME involves the propagation 
of three ground to space (uplink) laser beains as well as the reflection from space to 
ground (downlink) of one of these beams (see Figure 1). Two of the uplink lasers, one 
from AMOS (Air Force Maur Observatory Station). a site located at the Haleakala 
crater rim, and the other from Kihei. a site located at the volcano’s base, serve to align 
a relay mirror located on a payload orbiting over the Маш, Нажап stations. With the 
nurror aligned, à 1.06 u m laser travels from AMOS to payload and 15 relaved down to 
a hexagonal target board at Kihei.. The telescopes and optical detectors on the target 
board record the centroid motion and beam size of the reflected 1.06 u m laser beam 


Atmospheric turbulence has an important effect on both uplink and downlink Decis 


= 


(Nelson, 1988). Understanding and characterizing the atmospheric turbulence over the 


Hawaiian experiment site is crucial for assessing the success of the RME project. 


C. HAWAIIAN CLIMATOLOGY 

The following is a brief summary of Hawaiian climatology with a focus on the main 
experiment site in Kihei, Maui. 

The tropical island of Maui is located at 156? 15' longitude and 20° 45’ latitude 
(Blumenstock and Price, 1967). The island is the second most southern island in the 
Hawaiian chain and has an extreme length of 77 km, maximum width perpendicular to 
extreme length of 42 km and a maximum elevation of 3.1 km. Nicknamed "The Vallev 
Isle", this 1886 square km island is comprised of an isthmus linking two volcanos. The 
older volcano, Pui Kukui, is in the northwest sector reaching a peak of 1.8 km. To the 
southeast is the dormant volcano, Haleakala. Though Haleakala defines the island's 
highest point, it's crater is one of the largest on earth: 12 km long, 4 km wide and about 
900 m deep (Creal et al., 1987). Topographically induced high winds crossing the sum- 
nut are common. 

1. Climatic Regions of Hawai 

The tropical climate on Maui varies greatly with location. The dominate sub- 


climates include the following: 


NE un \IOUNTAIN slopes on the north, northeast WINDWARD SIDE (above 
600 m 


2. Moderately rainy and nuld WINDWARD LOWLANDS (below 600 m) to the 
merit, northeast; 
feet. LEEWARD LOWLANDS to the south. southeast: 


И ошегтлек raint, sometimes гу INTERIOR LOWLANDS between the two 
volcanic peaks; 


mae cold 1TIGIi NOL NTAINS above 600-900 m (Blumenstock and Price, 1967). 


Сл 


Rainfall measurements reveal the sharp weather contrasts over this island. On 
the leeward coastal areas such as Kihei, and mountain summits, average annual rainfall 
is 0.5 m. The lower windward slopes of a high mountain and at,near the summit of 
lower mountains (western Maui) have an average annual rainfall of 7.6 m (Blumenstock 
ШІСІ Price, 1907). 

2. Circulation 
Ihe dominant atmospheric circulation for the Ilawiiain islands is the north- 


northeasterlv trade winds. The prevalence of this flow reflects the seasonal north-south 
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movement of the North Pacific High. From May to September (summer). the North 
Pacific ieh is to the north. During this season, the southern sectors of the counter- 
clockwise flow around the High is such that trades prevail 80-95% of the time. From 
October to April (winter) the North Pacific High drops south such that the heart of the 
trades pass Well to the south of the islands. The frequency of trades during this season 
drops to 50-80% (Blumenstock and Price, 1967). 

Land and sea breeze circulation occurs on a much smaller scale. Kihei 1s unique 
in that land and sea breezes are present even during moderate to strong trades 
(Blumenstock and Price, 1967). The sometimes brisk sea breezes occur from forenoon 
ВЕЧНО evening. During the night. a very gentle land breeze flows until shortly after 
sunrise. The Kihei region also receives a nighttime Katabatic (downslope) flow compli- 
menting the land breeze alreadv in progress (personal observation). 

J. Temperature, Humidity, Cloud Cover 

Нажац“ mean annual temperature range generally fluctuates less that 6 degrees 
Celsius (Blumenstock and Price, 1967), Daily summer temperatures at 2.5 km are be- 
tween 4C (nighttime) and IOC (daytime). The winter range at this herght 1s 0-10C. In 
Memrowlands, daytime temperature is between 21-27C, while nighttime is in the 16-21C 
range. he varmest months in the dry areas such as Kihei are August and September 
(Blumenstock and Price. 1967). 

Generally humid conditions exist with high clouds except in the driest coastal 
regions and at high elevations. In the absence of svnoptic scale events. Kihei begins a 
day with locallv clear skvs. A backdrop of clouds clinging to the north and west slopes 
of THaleakalais common. Throughout the day Kihei experiences a growing and decaying 
of cloud cover. From sunset to approximately. 2000 hours local time. the partly cloudy 
Exeover over the site dissipates. Because of the dramatic climate changes over such a 
small island. clouds are generally visible at all hours from the Kihei experiment site in 
some section of the skv (Blumenstock and Price, 1967). 

4. Storms 

Generaly, lawanan weather falls into trade wind and non-trade wind condi- 
tions. During trade wind conditions one can expect trade wind showers on the E and 
NE (windward) coasts: possible light trade wind showers over the S and SW (leeward) 
areas; and, sea breezes ascending mountain slopes potentially producing light to occa- 
sionally heavy showers in the late afternoon and evening. Non-trade wind conditions 
are associated with light and variable winds; moist air slowly advecting onto the islands 


from the southeast, south or southwest; intense local storms; and, the possibility of 


~) 


major and minor storm svstems passing through the islands (Blumenstock and Price, 
1967) 


ОЈ 


Major storms for Hawaii fall into one of four classes of disturbances: 


Cold Front Storms - winter storm systems associated with cold fronts that ap- 
proach from the northwest end of Hawan, and move eastward barely reachmg 
Oahu. 


Kona Storms - winter storm svstems similar to cold front storms; winds are gener- 
ally from “kona” or leeward directions relative to the trades; well-developed kona 
storms are more widespread and more prolonged than the usual cold front storms. 


Tropical Storms - simular to hurricanes but with winds below 120 km hr. 


Hurricane - the least prevalent storm svstem: in 63 vears only four hurricanes have 
crossed over the islands. Unlike the Cold Front and Kona Storms, the Tropical 
Storm and Hurricane can occur at any ume. They are most likely to occur. how- 
ever, between the months of Julv to December (Blumenstock and Price, 1967). 


These events will be discussed in greater detail as they pertain to a specific case study. 


со 





ПІ. ATMOSPHERIC OPTICS THEORY 


A. ATMOSPHERIC INFLUENCE ON LIGHT PROPAGATION 

Light propagates through a medium in the form of a wavefront, "a surface over 
which an optical disturbance has a constant phase" (Hecht and Zajac, 1974). Fermat's 
Principle explams that the route light travels will be the smallest Optical Path Length 


(OPL). This is expressed mathematically as (Hecht and Zajac, 1974) 


Р 
ОРІ, -| n(s)ds, (1) 
S 


where n 1s the index of refraction (the "speed of an electromagnetic wave in vacuum to 
that in matter’), and s is the path length. The -phase of a propagating wave will be 
Кт даи times OPL (Eq. (1)). Optical turbulence arises when the accumulated phase 
along adjacent paths is different. Phase distortion in the atmosphere occurs when a 
Wavefront encounters random irregularities in the atmosphere’s index of refraction. 
Such irregularities are due to the turbulent mixing processes of the atmosphere. Thus, 
when a wavefront from a laser propagates vertically through the atmosphere, accumu- 
lated. random phase differences will degrade the laser beam and image system perform- 
ance. Ivpical atmospheric distortrons result in: laser beam centroid wander, 
scintillation, image breakup and blurring. 

The study of these atmospheric random turbulence processes requires the utilization 
of statistical quantities, means, variances and high moments. Frequently, variances and 
the higher order moments are not single-valued, unique or even bounded. Also, the at- 
mosphere is neither homogeneous nor isotropic... kolomogorov and Tatarski introduced 
the concept of structure functions and local homogeneity to resolve these problems 
(Tatarski. 1961). 

kohnogorov and Tatarski define the structure function. D,(r) , for an arbitrary 


spatially distributed random variable R(r) as 
DRE) = <IRE) - REDY >, (2) 


DUE -—r-—»,andtbhe angled brackets indicate an ensemble average of the enclosed 


quantity. Assuming the random process generating R is homogeneous and isotropic, 
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then only the magnitude of r is important, r = |7, — F,{.. Local homogeneity assumes that 
the mean <R(r)> = <R(r,)> (Chfford, 1978). Restricting the separation r to the inertial 
subrange of turbulence, Kolmogorov (1961) develops a universal form of the single 


Structure function 
Ре) я Сік, (3) 


where C} isa structure parameter, a measure of “the total amount of energy in the tur- 
bulence” (Clifford. 1978); ris the distance between the two sampled points; and p ime 
be between the turbulence innerscale, 4. and the turbulence ошегзсаје, 1, . Еда на 
below the innerscale correspond to the eddv's energv dissipation through viscous effects. 
The outerscale corresponds to "the largest scale size for which the eddies mav be con- 
sidered. to be isotropic (СПП vo 

For optical propagation through turbulence, the characterizing quantity is the index 
of refraction structure parameter (Walters and Kunkel, 1981) 
y. “(= mc 

s — F ' (4) 
p 

where a = n(rj) and лп, 2 n(r,) are the atmospheric index of refraction at points K 
respectively. Though C? is difficult to measure directly, profiles calculated mine 
rawinsonde sensor measurements show that C? generally decreases with altitude (Nelson, 
1988). A more pracucal approach for measuring atmospheric turbulence 1s to measure 


the temperature structure parameter С3. 


2 «cm E ше Я 
Ст = sa s. 0) 
y 


Since the atmospheric index of refraction 
EUR (6) 
Г 
depends on pressure P and temperature T (Tatarski, 1961), C? 1s related to C5 bv sub- 
ѕпоппе Ед. (Оо Еа оза 


ге дии 


^ 


= (79 х 10" ==) Ст. (7) 
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assunung turbulence 15 150Багіс. Specific methods for calculaung C2 will be discussed in 


Chapter 5. 


B. OPTICAL VARIABLES 
As described earlier, the laser beam 15 affected by atmospheric turbulence as a result 
of index of refraction fluctuations along the turbulent propagation path. Two parame- 
Memes that characterize the optical turbulence generated by the atmosphere are: the 
Bunsverse coherence length, 7,, and the isoplanatic angle, 6, . 
1. Transverse Coherence Length 
ИО рос length 7an sa measure ol the lateral autocorrelation 
length of the electric field propagating through the turbulent atmosphere. Because of 
the Fourier transform properties of an unaging system, и, 1s also a measure of the spatial 
frequency response of the atmosphere. the atmospheric modulation transfer function 
(MTF). Viewed as a component of a linear system, the atmosphere acts like a low pass 
filter suppressing the high spatial frequencies. or details of an image. 
The perturbation theory developed Бу Tatarski and extended Бу Fried (1966) 
and Lutomursk: and Yura (1971) shows that the long term. modulation transfer function 
of the atmosphere is expressed in terms of the wave structure function D bv 


NN р) 
MTF(v) = exp( —— ), (8) 


where v 1s the spatial frequency in cycles m. 4 is the wavelength of light and fis the focal 


length of the opucal svstem. For a plane wave, D has the form 


ја 
Cia. (9) 


D(p) = 2.911202 


EU 


where k is 27/7, the wavenumber; p = |F, — r,| + fv, a distance measured at the entrance 
pupil of the optical system. L is the optical path length and C is the index of refraction 
structure parameter. 

The general format of MTT(v) 15 exp( —ap**) where a is a turbulence and wave- 
length scaling parameter. MTT(v) is frequently written as: 


- 


ІЙ en | 
МТД) = exp| — 3.44 = сехр = 344 P Pi (10) 


representing Fried's (1966) version of the long term atmospheric modulation transfer 
function. Equivalently, Lutomirski and Yura’s (1971) mutual coherence function (MCF) 


15 
МСКр) а екрі зо У | (11) 


where p, 1s an autocorrclation length of the electric field. 
The parameters x, and p, are different measures of the electromagnetic (EM) 


wave s coherence as restricted bv the atmospheric turbulence. Thev are related by 


» 


Using Eq. (8) and (9). the coherence length for a plane wave, ry, becomes 


L 
го = 2.1[1.46k? | сауар?” (13) 
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The parameter ry represents the distance traverse to the direction of propagation where, 
on average, the correlation of the electrie field ise o ЕТ” 
To interpret ry consider the following: 


l. For an optical aperture d smalier than the coherence length. the electric field across 
the aperture is colierent.. Tlie angular resolution will be proportional to 4/4. 


tJ 


l or an opucal aperture larger than the coherence length. the correlation across the 
aperture is limited to regions of ry or py (depending on how tlie correlation 1s de- 
fined). fhe LM intensities across the aperture ушШ ач (епегех 15 conserva l 
opucal svstem's average angular resolution will be proportional to 2o C 
Favier and ОНА 


(о 


Visually. for apertures smaller than 5. atmospheric turbulence results in image 
centroid motion, a process called, “tlt”. In telescopes larger than 7, the image 15 
broken up into muluple speckels. each associated with a 7 sized region. The angle 
A/r, determines the overall envelope of the complex image. 


4. With respect to measuring the atmospheric optical turbulence, large г, values, such 
as 200-400 mm, indicate a small amount of optical turbulence is present along the 
integrated path. Conversely, small x, magnitudes, such as 20-40 mm. unply a large 
amount of optcal turbuleice exists along the integrated path. 

2. зорапайс Апо!е 
lhe isoplanatic angle, 0, , is an angular measurement of spatial coherence be- 


tween two intersecting ravs. (Walters, 1983). Given two liglit sources a small angleag 








apart, the random atmospheric phase fluctuations along the two paths are different. 
Figure 2 displays ravs r,, and r,, interacting with three different turbulence patches. At 
the atmosphere top the two paths undergo unique path distortions. At the surface (near 
point 1), the turbulent interaction is common to both rays. 6,15 the angle between rays 
га and r;, such that the accumulated phase fluctuations are correlated to within e! of a 
perfect correlation. 


Fried (1982) develops a mathematical representation for @, for apertures d > > 


Ға. 


ОФ 


бо = [2.988 h s 31, (14) 


ШОО 22/2 wavenumber of light; C? is the refractive index structure parameter; and, 
2 15 the alutude from the ground. 

beca c osolles displav a rapid decrease with altitude. The 255 terin is sig- 
nificant in that it emphasizes the optical turbulence at higher altitudes. The net product 
of these two terms (the integrand) results in a profile of relatively constant values from 
mew kilometers to 10 — 15 kilometers (Walters. 1965). 

Wie practical interpretation of @, 1s simular to 7,. When @ magnitudes are small 
(1-2 u rad). the optical turbulence along the integrated path and especially around the 


tropopause level is large. Conversely. when 0, values are large (14-20 џ rad). optical 


— 
` 


turbulence 1s small. 


Ро. 2. Diagram of Isoplanatic А Ис 








ту. ВМРЕ МІ В ТАРУ ДР А 


A. RME KIHEI DATA COLLECTION 

The initial sampling of isoplanatic angles, 0, took place on the summit of Mt. 
Haleakala in Maui, Hawai. These measurements commenced on 17 September 1985 
and concluded on 1 February 1986. Simultaneous meteorological surface measurements 
at the site provided some understanding of localized atmospheric events. A vear and a 
half later, both 0, and r, were measured at the base of Haleakala (at approximately sea 
level), in Kihei. Six sampling sessions of four and five mights took place between August 
1957 and July 1988S. Two Hawai Rawinsonde Observation stations, [lilo and Lihue, 
provided raw data hstings for each of the 1987-1988 cases. Additional on-site 
rawinsonde launches supplemented the optical data for four out of the six sessions (Ta- 
ble 1 ). The meteorological data for the August and December 1987 experiments used 
the VIZ WO-SOOORP+ rawinsonde svstem with 10-second vertical resolution. The 
April and July 1988 experiments emploved the Vaisala Digi CORA MW 11 system with 
2- and 10-second vertical resolution. All four cases used the OMEGA Navigational 
Network for wind measurements. Intermittent GOLS-WEST Infrared Satellite Images 
complimented the rawinsonde data for all the cases. The primary focus of this study 
will be on the April and Julv 1988 experiments, due to the higher vertical resolution 


provided by the Digi CORA dataset. 


Table 1. MAIN OPTICAL AND METEOROLOGICAL RME SENSORS 
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B. KIHEI RME EQUIPMENT 
1. Optieal Equipment 

Optical measurements were taken by an rsoplanometer and transverse coherence 
length sensor desrened burlt by Dr. D.L. Walters, NPS. For spectfre details on састы 
strument the reader should pursue Stevens (1985) for the isoplanometer and Walters, 
Favrer, and Hines (1979) for the transverse coherence length sensor. Isoplanatic angles 
were measured bv an isoplanometer mounted on a portable 200 mm Celestron telescope. 
The software was run on an HP 217 computer. Isoplanometer data were automatically 
sampled once every second and displayed real-time on a time verses isoplanatic angle 
plot. Before storing the data, the samples were converted mto 10 ѕесола ахегасез а 
complete isoplanatic angle system 1s capable of operating independent of human inter- 
action for 1-3 hours provided polar alignment is perfect and the chosen star's zenith 
angle does not exceed 40 degrees. 

Тре transverse coherence length sensor was mounted onto a portable 355 mm 
Celestron telescope. Software used for real-time processing required an operator to 
evaluate the image and stellar positioning before each sample. This added interface with 
a trained researcher served to suppress instrumental artifacts. The average sampa 
was between 1-2 minutes sample. Processed data were stored on anillo P E: 
computer. 

2. Meteorological Equipment 

Meteorological measurements specific to this thesis were a product of the 
Vaisala DigiCORA MW 11. This solid state. digital design rawinsonde package mte- 
grates the receiver, data processors and operator display into one unit. High speed 
precision sensors pernut rapid sampling (at approximately 2 second intervals) resulung 
m high density. data acquisition. The 95 mm by 145 mm (approximate size) RS 80 
Rawimsonde (Navaid Radiosonde) package was launched on an unravelling tether at- 
tached to a 200 gm balloon filled with helium. Pressure values were measured by a 
BAROCAP pressure sensor, a steel capacitive aneroid with a resolution of 0.1 hPa and 
an accuracy of 0.5 hPa. Temperature sampling was taken by a [HIERNIOCAT ола 
a small capacitive bead in a glass casing. Resolution of the TIIERMOCAP was Q.1C: 
accuracy was 0.2C.. The HUMICAD, a thin film. humidity sensor. measured dmm 
with a resolution of 1? relative hunudity (RII) and an accuracy of 276 RII. Accordmg 
to the Vaisala Ref. A0435. this sensor produced a reliable response even mn low temper- 
atures and after exposure to condensation. The measuring range published was 0-100" 
IF (Vaisala; 19563): 
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The commutation between the sensors and the transducer electronics was ac- 
complished by a digital electronic switch. А water activated battery provided the RS 80 
“sonde” with a light weight power supplv capable of keeping the package activated for 
the entire surface to 30 km (balloon permitting) launch sequence (Vaisala, 1986b). 

The telemetry frequency was between 400-406 MHz. Wind data for this exper- 
iment were calculated from the OMEGA Navigational Svstem. The RS 80 received 
signals from the OMEGA Navaid Network (11-14 kHz). superimposed this information 
onto the thermodynamic data and relayed the sampling to the DigiCORA MW 11 
(Vaisala, 1986b). 


C. KIHEI SONDE DATA PREPARATION FOR ANALYSIS 

The initial form of the DigiCORA MW 11 rawinsonde output used in this research 
began as a hardcopy listing of time in minutes and seconds, height in meters, pressure 
in millibars, temperature in Celsius, relative humidity in percent, dewpoint temperature 
in Celsius, mixing ratio in grains per kilogram, wind speed in knots and wind direction 
EMUeerees. Using а Dest PCSCAN scanner. tie powerful IBM mainframe XEDIT. a 
software package capable of transfering IBM files to a IBM PC-AT, and the instant 
graphics displav of QUATTRO, the initial hardcopy was duplicated as an ASCII file. 

Thermodvnanuc sampling rate for the first two kilometers was every two seconds. 
Bevond this threshold. sampling continued once every 10 seconds. Due to the limita- 
tions in the Omega Navigational system, winds were measured every 10 seconds 
throughout the launch. A hnear interpolation between the 10-second wind measure- 
ments was applied to the initial two kilometers remedving the inconsistent sampling rate. 

Variables calculated from the raw data include potential temperature. potential 


Exuperature eradient, shear. Richardson number, C2, C2. n; and 6,. 
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V. ANALYSIS/RESULTS 


lo strengthen the link between the measured optical parameters апа Ше 
meteorological data collected. a number of methods for calculating r, and @, values were 
pursued. As explained earlier, both r, and 8, are a function of C? , which, in turn, 15 а 
function of C? (see Chapter 3). In this chapter, four methods for calculating C? from 
rawinsonde data will be investigated. However, the first step is to define an appropriate 


conservative temperature quantity to be used in the C? calculation. 


A. STANDARDIZING THE TEMPERATURE PROFILE 

An air parcel vertically displaced in the atmosphere experiences a temperature 
change as the parcel’s pressure equalizes with the ambient environment (Clifford, 1978). 
Permitting the air to be moist but unsaturated, and assuming the atmosphere 1s in 
hydrostatic equilibrium, the rate at which the lifted parcel’s temperature will decrease 15 
the drv-adiabatic lapse rate. This 1s approximately 9.767C km. If the parcelis satum aa 
the release of latent heat produces a lapse rate that is less than the drv-adiabatic lapse 
rate (Iiuschke. 1970). For this study. it is important to work with conservative variables, 
a variable that is independent of position in space. As described above, temperature 
does not fit this description. The conservative quantity (independent of position) best 
suited is potential temperature (Clifford, 1978). By definition, potential temperature 15 
“the temperature Which the parcel of air would have if. it were expanded or compressed 
adiabatically from its existing pressure and temperature to a standard pressure p, (gen- 
erally taken as 1000mb»' (Wallace and Hobbs., 1977). In short. the temperatures sam- 
pled vertically are standardized to one pressure level. Mathematicallv, the potenual 


temperature expression is 
оа | 
©) — ТІ ) С, ° (15) 


where © is potential temperature in kelvin. T is temperature in kelvin. /’y is a standard 
pressure such as 1000mb, P is pressure m millibars, R is the gas constant for dry air and 
C, 1s the specific heat of drv air at constant pressure (1.005 Jg K) (Huschke. 1970). 


Note: Following the usual convention, Cj shali hereafter be written іш ріасе оГ % 
(Chflord, 1978). 
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ER CALCULATING C}: FOUR CASE STUDIES 
1. Case E: Simple Potential Temperature Differences 
Case I takes the given mathematical definition, 


ee ОЗ О) = (О ог» (16) 
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and assumes ©,=0©,. Physically, this presumes homogeneity between neighboring 
samples. In other words, the samples are taken within the inner and outerscale of the 
Eu ent layer(latarski, 1961). Note: the angled brackets indicate an ensemble average 
of the enclosed quantitv. Associated with the above and all the cases hereafter is the 
assumption that the sampled laver 1s isotropic (Tatarski, 1961). 

Unlike most C} calculations, the dataset employed by this research consists of 
vertical temperature samples. For the above assumption to be true, the rate of sampling 
must be sufficiently high with respect to the potential temperature gradients. If the 
sampling rate 15 too slow, trends in potential temperature will show up in the differences, 
(O,—©,). Figure > displays such trends first between 2 and 5 km then again above 12 
kin. Note: The large increase in potential temperature differences in the stratosphere 
was found in all the profiles examined. 

a. Results 

sine ihe Cite potential temperature profile for the C? calculations, the 
eight April 1958-Kihei rawinsonde flights and the ten Julv 1988- Kihei launches all under 
estimate the direct r, and 0, measurements bv a factor of ten. (See Appendix A) This 
15 expected considering the trends in the potential temperature differences. As demon- 
Strated in Figure 3, these trends generate differences approaching 2 K near the boundary 
u In the stratosphere (above 15 km). potential temperature differences are even 
larger. For this particular profile, the stratospheric potential temperature differences 
reach as much as J K. When these large differences are blindly used to calculate C}. the 
result is a deceptively large over estimation of Cz. The over estimated C$ subsequently 
leads to an under estimation of r; and 0. 
2. Case I]: Using a Running Mean 

Case I] investigates the effects of detrending the potential temperature profile 

prior to calculating Cz. The two representative high pass filters selected for this case 


include the equally weighted 11] Running Mean (RM) 


ISS: 


{ 


IA а а, 


ә 


(17) 


and the weighted 13531 Running Mean 
у= Ір ок ЗА OA не и а | (18) 
13 
Figure 4 displavs the potential temperature difference profile after the 111 RM has been 
subtracted from the original potential temperature profile. There is still an increase in 
potential temperature differences around the boundary layer, however, the peak values 
of approximately 0.3 K is balanced with an adjacent value of about -0.5 K. Even the 
stratosphere appears nearly Svminietnical betwiecrmaamen E 
а. 111 Running Mean 
The implementation of the equally weighted (111) RM handles the 
DigiCORA's two- to ten-second sampling rate transition bv using values equally spaced 
in time when available. When data are unavailable. preference is given to tlie Commas 
point. Profile endpornts are ignored. For two-second data the filter takes into account 
layers between 15-30 meters thick. Ten-second data involves Jlavers of 90290 mem 
b. 13531 Running Mean 
The 1383] running mean filter accommodates the two- to ten-second sam- 
pling rate change with preference given first to data equally spaced in time. and second, 
to the central point. Endpoints are noted, but generallv ignored. 
З 
Both RM methods produce к values within the same order magnitude as 
the optical sensor values. With few excepuons, the relatively large and small magnitude 
variations present in the optical sensor data are traced bx both RM гезше. Лени 
RM generate values approximately 16? larger than the 111 RM. The error between the 
111 RM and sensor results ranges from 13% то 55%. ТПһе 13531 КМ опас еи ни 
the sensor values by 4°0 to 48%. When the C profile detrended by the 111 RM is 
truncated just below the stratosphere's large C? increase, r, increases bv only 11%. "fae 
ble 2 displavs a sample of these results. 
d. 0 , Results 
While large scale magnitude changes recorded by the sensor are virtually 
ignored, the calculated @, values are within the same order of magnitude as the sensom 


m 


The range of sensor measurements 1s from 3 to 7.5 u rad. Values caleulated (топ та 
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111 КМ cover a range of 0.7-1.3 и тай. (Те 8.5 yw rad of flight 5, July 88 is ignored.) 
EealculatcdeviththeslsssIeRM rangesefrome0.8-1.5 grad. Truncating the detrended 
C? profile just below the stratosphere's large C3 increase, increases 0) by 7820. A sample 


of these results 1s in Table 2. 


Та е 2. SINGLE PROBE RESULTS - 0605 UTC 28 APRIL 19883 
ПРАГ СОВА ЛОХ 


SPTD (all values Eu — 
Кене | _ 


BULL I VE (all values) 


111 RM 
{troposphere only) 


1353] RM (all values) 
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3. Case III: Using the Bulk Gradient C; Equation 
The third approach incorporates an expression devcloped bv Tatarski (1961) 
Euemlates C5 to bulk parameters of the atmosphere. This expression which will here- 


Ире referenced as the "Bulk Gradient C? " equation (BG) is. 
Gu (O jj). ШШ) 


where a is a constant: х 15 the ratio of eddv thermal diffusivity to eddy viscosity; L, is the 
outerscale length: and. (dO/dz) is the potential temperature gradient. Because L, 1s not 
universally well-defined and dO/dz is also vague with respect to the resolution needed to 
property define the derivative, the values calculated from BG are suspect. The main 
Strength of the expression is that it accentuates the strong dependence C} has on the 
square of the temperature gradient. This term also provides the "atmospheric density 
inhomogencity susceptible to being mixed bv the wind turbulence" (Walters and Kunkcl. 
1981). Assumptions made when applving this equation include the following: 


1. The atmosphere is mechanically turbulent. Note the lack of a shear term. 


гә 


The outerscale, Z4, is half the distance between samples 
Шарл, 215. 2.0. 


ЖЕК апо. ж 15.1 . 
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а. Results 

The BG r, and @, results are one-tenth the sensor measurements. Both op- 
tical parameters calculated in thus case closely trace those values generated by the Simple 
Potential Temperature Differences (SPTD) Case. r, magnitudes vary over a range of 6-8 
mm. @, values cover a range of 2-3 mm. These values are within 7% of the ЭРТІР СЕ 
Considering the dominant term of the BG equation is the potential temperature gradient, 
it is not surprising to discover the finał results are so closely associated to Case I. An 
explanation for why r, and 0, Values are off by a factor of 10 begins with a re-evaluation 
of the assumptions. The constant, a, has a suggested range between 2.8 to 3.2. Shifting 
this to the higher constant will not accommodate the discrepancy. Changing the 
outerscale length to the full thickness of the laver serves only to decrease и папи 
Justifving an outerscale smaller than half the thickness of the laver sampled taps into the 
fundamental problem with the C} calculations. That is, vertical resolution. The sampled 
values are too far apart to provide the temperature gradient information needed to drive 
this C2? eyuation properly. Decreasing the outerscale length without simultaneously de- 
creasing the thickness of the sampled laver 1s illogical. As in Case I, the BG results in- 
dicate a need for either detrending the data prior to calculating C} or using a rawimsonde 
dataset capable of producing more detailed temperature gradient information. 

4. Case IV: Using Richardson Number 

Сз. 15 a constant representing the total amount of energy in a turbulenta а 
(Chfford. 1978). The vertical temperature samples measured bv the rawinsonde авевиии 
always taken from turbulent lavers. To computationally isolate onlv those atmospheric 
lavers which are turbulent, the Richardson number. Ri. was used. Richardson niuimber 
is “the ratio of Work done against gravitational stability to energy transferred from mean 
to turbulent motion” (liuschke, 1970). When Rr < 0.25, the laver 1s unstable ande 
calculated using the SPTD. Running Mean and BG approaches, respectively. When Ri 
> 2. Ше laver is stable and no C} calculation is made. “Critical conditons Crio 
0.25 < Ri< 2.0. As the final case reviewed, Cis calcule dEn e ЕК ни 

а. Results: Unstable Ri Case (Ri < 0.25) 

Letting a turbulent laver be defined as one having a. Ri « 0.25, each Olam и 
above three methods for calculating C was conditionally computed. With remarkable 
consistency, all three methods over estimate r,.. y values farthest from the sensor's value 
generally are produced by the SPTD case. Except for one eXtreme сазе леге ты 
culated as greater than 300 mm from the sensor value, the SPID r, хае аге 


2 to 3 times too large. The BG method leads to values that are approximately 7° closer 








than SPTD to the sensor value. The RM method(s) r, values over estimate the sensor 
nasto 71mm forthe 1353] RM and 8 to 58 mm for the 111 RM. 

0, calculations generate less consistent results. The sequence of greater to 
smaller magnitudes seen in the r, results does not hold for 0,. About 40?» of all the 0, 
Eu calculated for Case IV (R1« 0.25) over esumate the sensor's measurement: 5095 
under estimate б„; and, 1070 straddle the sensor value (see Appendix A). 

b. Results: Conditionally Unstable Ri Case {Ri < 2.0) 

With the inclusion of more atmospheric lavers in the r, integral, one would 
Sect tne resulting magnitude to be smaller. For the most part this is truc. The SPTD 
Ehe BG cases respond in the same manner. Their 7, values under estimate the sensor 
ШІ 59% Гог Ше ӘРТІ? сас апа 22-65% for the BG case. The 111 and 13531 RM 
calculations still over estimate rj, with the 13531. RM slightly higher than the 111 RM 
(See Appendix A). The magnitude of the results as compared to the sensor's value is 
Benerallv within 7826 for 13531 RM and 59% for the 111 RM. 

@, results display a pattern consistent with the “All Values” results. The 
u 10 Smallest 0, magnitudes calculated are: (Sensor Valuc), 13531 RM, 111 RM. 
отапа 80. АП хаџев аге 1255 than 5 u rad! As before the large magnitude vari- 


ations observed in the sensor data are virtually ignored in the calculated data. 


С. CORRELATION OF SENSOR AND CALCULATED PARAMETERS 

Correlating the measured and calculated optical parameters for cach case discussed 
КИ ООП D senerates r-squared values between 0.0009 and 0,5. Table 3 lists r-squared 
ШЕЕ Гог each case. Note: r-squared is | for a perfect correlation. 

See lack Of соод correlation demands.a review of the sensor's ability to sample 
values for C3 calculations. After detrending the potential temperature profile with a 111 
emer potential temperature’ differences for the two-second data are plotted. (See 
Figure 3) Initially, this profile appears to contain valid turbulence information. Re- 
peating this procedure with the temperature (K) profile (Figure 6). however, unveils a 
major problem. The systematic + -0.03 K oscillation dominating the profile implies 
that the temperature differences below 2 km are predominantly quantization noise from 
the 0.1°C temperature resolution. Comparing Figure 5 and Figure 6shows that the 
potential temperature calculation introduces additional errors that obscure the quantiz- 
ing error. Rather than pursue a detailed correlation analysis. attention is diverted to 1) 
identifving and quantizing the major sources of error, and 2) providing recommendations 


for a temperature sensor system more suited for calculating C2. 
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Fig. 5. Initial 2 km of Detrended Potential Temperature Dilferences 
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Table 3. MEASURED AND CALCULATED PARAMETERS CORRELATED 


R-SQUARE ! 
ee ro (mm) Ө, (pg rad) 


Ce CALC 


LATION ` a 
METHOD | 


ДЕ 


x 
Ri «0.25 Values 





D. ERROR ANALYSIS 
l. Temperature Probe Time Constant Correction 
In both the April SS and July 88 soundings, the temperature was measured with 
a capacitive bead. For sensor specifications see Chapter 4. To calculate the sensor ге- 


sponse time c , the following equation dss 
ре ее ! m ~ 
ша i Al a a D. b * (20) 


where 7 15 velocitv. A is the constant for the calm atmosphere (v 2 0) and B 1s a constant 
at velocity. v. For this particular probe, the sensor time constant at a velocity of 6m 5 
15 2.5 s. With an ascent rate of 4 ms, t increases to approximately 3 Ss.” Ш ой 
calculate the error induced bv the temperature sensor time constant, the structure func- 
tuon of the probe system with limiting scale lengths must be compared to the structure 


опсом охота Ссс 


Mathematically, the structure function over all scale lengths 15 (Tatarski, 1961): 


D(r.0,00) = 4 | sin’ (kr/2)k dk, (21) 


ihe 


г 


where k is ? z 4.27 the turbulent wavelength and r the probe separation distance. 


The structure function filtered by a one pole, low pass filter is (Hunten, 1964): 





b 
sin^(kr] 2)k 77? 
Драйв) 2 4| ——,——ak, (22) 
| + (от) 


кеге 
а 15 2z/ L4, the limiting lower frequency (outerscale, L,), 
bis 2z/l,, the limiting upper frequency (inner scale, 4). 
с 15 KV where V is the wind velocity (rawinsonde’s ascent rate) and 


t is the time constant of the temperature probe (the upper frequency cutoff). 


Тһе от term in the denonunator attenuates the high spatial frequencies. Since the k-*? 
Eununcreases at low wavenumber, Eg. (22) 1s much more sensitive to the selection of 
an outerscale than the innerscale. Test calculations of the D(Sm.a,b) integral confirm 
that an innerscale of 0.1 m verses 0.001 m produce virtually identical results. When 
varving the outerscale from 10 m to 1000 m (all other parameters of D(Sm.a,b) remain- 
ing unchanged). the 1000 m structure function is 200 times larger than for the 10 m. 

With an inner scale of 0.1 m. an outerscale of 1000 m and a 3-second temper- 
горе time constant, the DigiCORA’s 2-second data 1s 1596 of the actual C 
present. The structure parameter C3, calculated from the 10-second data is 54% of the 
actual C? value. While the 10-second data appears to produce results more than three 
times closer to the desired values, this is somewhat deceptive. The larger time interval 
between samples (accent rate is constant) results in a greater probe separation distance. 
Consequently. higher spatial frequencies or rapid changes in the potential temperature 
will be missed. This aliasing of temperature measurements 1s one of the major difficulties 
with calculating structure function parameters from rawinsonde data. 

Another problem with calculaung a structure function D from 10 second (r= 40 
m) data involves a limitation of the Kolmogorov r° slope model. D is a function of the 
Probe Separation distance, r, and the outerscale. For r values comparable to the 
outerscale. the actual D(r) deviates from Kolmogorov’s model. When r is larger than 
the outerscale, the measured (actual) structure function will be smaller than the D cal- 
culated from Kolmogorov’s model. (See Figure 7) The significance of this problem is 
that the calculated D under estimates the actual turbulence present when r is greater 
ieetne outerscale. Thus. accepting the calculated C2 as 54% of the actual C} present 


15 done cautiously. 
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Fig. 7. Kolmogorov Model vs. Actual Structure Function Curve. 
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а. Application of the Filter 

Correcting the basic C7 equations for the temperature probe time constant 
produces an even greater over estimation of atmospheric turbulence. The magnitudes 
of the filtered қ апа 6, (all values) decrease 35-50%% from the unfiltered values. 
Table 4 and Table 5 display the results taken from two rawinsonde datasets: 0605 UTC 
28 April 1988 and 1434 UTC 30 April 1988. Imposing this correction onto the calcu- 
lation of the optical parameters where R1« 0.25 and Ri<2.0 produce the same 35-50% 
decrease as the "All Values" cases. Unless otherwise noted, further reference of cor- 
rection implementauon will be done on the 0605 UTC 28 April 1988 Kihei, Маш 


rawinsonde dataset. 
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2. Sensor/Atmosphere Interaction Noise 
Since the optical parameters calculated bv the rawinsonde are smaller than the 
optical sensor values and compensation for the probe response-time accentuates this 
problem. some other large error source must be present. Noise introduced by the 
sensor-atmosphere interaction 1s such a source. Two methods attempt to quantize this 
phenomena. 
а. Standard Error Formula 
The first approach tackles the problem bv focusing on the limitations of the 
rawinsonde svstem. [n other words, knowing the resolution of each variable measured, 
how well can the given instruments measure atmospheric turbulence? This method uses 
the standard error formula which for a function E(x,v.z) 1s (Bevington, 1969): 
C Tono ie CE MEET 3 


СЕ = (=) (an nei m ) (у) 2 Фе rear (23) 





Applying the standard error formula to the SPTD calculation of C? produces 


ТИ E 2 чт Ty DNE x 
= Ј он р: а cssc í yn с. l. (24) 
E лк 


Note: 7, and 7, are assumed to generate independent random errors. 
With C denned 


C9 x10" = аа (25) 
the svsteml induced eiemfor C715 
= e 2 š 
- ОРОШ A 7 асое пара 
Дан = Е | (Pacz) ES (т) ЕЕ (2Суср) . (20) 


Implementing the error correction reveals three layers where (C7 — боз) <Q. Theseus 
sults imply that the measurement uncertainties are comparable to the measurements. 
The sensitive regions include 0 to 2 km, 4.5 to 9.5 km and 12 го 14 кт. s и 
area of major concern is near the surface where the majority of turbulence 1s measured. 
Eliminating the 0-2 km laver would result in a gross under estimation of turbulence 
which ultimately generates large sy and @, . Since by definitron C? can not be less tham 


zero, all negative error corrections are assunied TosDerzenas 


о 
tu 





The results from the C? error calculation show values less than or equal to 
lxccarthe surface (0 - 2 km). Between 3.5 - 9.5 km, aswell as. 10 - 14 km, the C? < 0 
condition 1s intermittent. Unlike the C3 calculation, C2 has a pressure term that effec- 
tivelv minimizes the influence of any values above the tropopause. Thus, the main 
concern is with turbulence measurements near the surface. 

Пе < пистите parameter, including sensor error (0605 СТС 28 
295 1955), 1$ (1.1 - 1.8)х1Ю-3А?м-23 Гог the 0 - 2 km two-second data; the average 
ЕЕ Сеп5от еггот отег this same region is ($ -- 13)x10-7155:15 , The 10-second data (2 km 
and above) produces an average C+ + sensor error of (7.4 + 0.8)x10-7A7»n-?; the average 
Маса по sensor error for the 10-second data, 1s (2.0 4- 0.4)x10-5pr??. Clearly, the 
0-2 km data are buried in uncertainties. 

b. Vertical C} Quantization Noise Simulation 

The second technique used for quantizing the roundoff error in the SPTD 
(Simple Potential Temperature Diflerences) calculation of C2 simulates the atmosphere's 
temperature profile with a straight line slope of -0.007 K'm. This technique (Walters, 
1989) includes a correlation between adjacent temperature samples that 1s ignored in the 
standard error approach. The 0.1 К temperature resolution produces a quantization 


D 2 


noise of I.-1x1024A7m^?? for the two-second data and б6.5к10-А?т-? 3 Гог the ten-second 
data. The two-second noise 1s 22%o less than the Standard Error Function's result. The 
ten-second value is the same. Running the simulation on a detrended potential temper- 
ature profile. the calculated noise for C3 is a factor of 10 smaller than the SPTD calcu- 
lation. Specifically, the two- and ten-second quanuzation noises are 6.2x10-^A?5r?? and 
emery Anc? *, respectively. 
3. Implementing Noise and Filter Corrections 
ite measured C+ is 


е: е (СТ > Ст) + Noise, (27) 


measured 


Where the Filter is the temperature probe time constant correction, and Noise is the 
vertical C} quantization noise from the srmulation. 

Applying the noise correction to the 111. RM C3 calculation for the 0605 UTC 
28 April 1988 rawinsonde flight, the original 36 mm 7, calculated by this method in- 
creases to 46 mmm (all values included). The sensor value is 48 mm. Replacing all neg- 


ative C3? with zero reduces this measurement to 41 mm. Imposing the temperature probe 


за 


time constant correction reduces this value again to 36 mm (all values included) and ul- 
umatelv 7, 15 25 mn including onlv positive C? values. 

0, displavs a much smaller range of resuits due to the equation's sensitivitv for 
values at the tropopause coupled with a lack of large variation created bv the noise and 
filter at that level. Specifically. 0, starts as 1.3 и rad (only the 111 RM applied). In- 
cluding all values and subtracting out the quantization noise increases 0, to 1.4 д гай. 
Removing all negative C} does not change the calculated measurement. Dividing the 
noise-free C+ by the filter leads to a value of 1.0 и rad (all values). And finally, using 


only positive C3, the 6, result remains 1.0 un rad. The sensor's measurement 1s 3 д таа. 


E. DETERMINING AN ADEQUATE TEMPERATURE RESOLUTION 

The (0.1 + 0.05)K temperature probe resolution dominates the error termo ин 
(24). An mportant consideration is to determine the temperature resolution required for 
calculating C$ from rawinsonde data. First, let oe, = 1xIO-*&m??* . А C3 noise level of 
1x 1072 Ar? is adequate for most measurements. Assuming the temperature difference 
is 0.0098 K over a 1 meter height difference (drv-adiabatic lapse тасс). ачи al 
с., + 6, = 0, опе сап solve for c}. Note: 05, = 07. The result is c} = 1.30x 10* A9? 9. 


The roundoff error (noise) temperature variance expression (Bose, 19835) 
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2 -24 

сү О (25) 
where dis the number of base ten significant digits to the right of the decimal point re- 
maining after roundoff, permits the calculation of significant digits required for a given 
variance. Solving for d. the resolution necessary for a С съ. == Іх а арргохі- 
mately two десила! роии$ гетаииие аКег гоипЧойЙ, 0.01К гезомшиоп. Repeating the 
above procedure with a А Г/А: of -0.007K 1m (a lapse rate observed in tlie Kihei data- 
set), the results are also approximately two. 

The initial use of the standard error formula defined o, as 0.05 Kk, as implied bv the 
sensor specifications. The roundoff error temperature variance expression calculates o; 
as 0.03 К. When 0.03 К is used in the C} standard error formula, the product 15 ap- 
proximately one-third the magnitude of the standard error formula using 0.05 Kk. This 
decrease in the C$ noise term increases C}. A larger C} leads to a larger C? and, Subse= 
quently a smaller 7, and @,. Since the calculated optical parameters are already under 
estimating the sensor values, this adjustment only confirms that the current temperature 


sensor svstem 1s а: best marginally acceptable for calculating C;. 





F. RUNNING MEAN FITIING ERROR 

A fitting error around sharp gradients, such as across inversions, arises when im- 
plementing the running mean into the €? calculation. The running mean is unable to 
follow the trends. Figure $ displays the large difference that occurs between the running 
mean and a given profile. Though variances are used in the C2calculation, the variances 
@revartificially enlarged Бу the poor fit across the inversion. The SPTD experiences a 


form of this fitting problem that results in aliasing. (See Section DI of this chapter). 
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VI. DISCUSSION AND CONCLUSIONS 


A. CORRELATED OPTICAL DATA AND METEOROLOGICAL EVENTS 

Optical and meteorological data for the Relay Mirror Experiment (RME) were col- 
lected between 1985-1988 in. Kihei, Maui. Hawaii. Reviewing the composite Kihei op- 
tical dataset reveals the Maui site as “mediocre” to “poor”. While upward and downward 
trends occur, the predonunant transverse coherence length values generally stay within 
Neto /Ü mm range. 

The lack of extremelv good values (e.g., 200-400 mm) is tied to the turbulent, 3-5* 
C boundary laver inversion 2-3 km above the Kihei site. Except in the winter, north- 
northeasterly trades usually cap the boundary laver. The interface between the 4 to 7 
m.s trades and the boundary laver provides the mechanical turbulence detrimental to 
optimum “seeing” conditions for hght propagation. Ihe concept of the boundary laver 
introducing optical turbulence is supported by Gossard (1985). Using fast-response 
turbulence sensors traversing a 300m tower, he shows a turbulent laver of approximately 
ENSE? to exist just below the inversion. This implies that the greatest mixing of air 
masses with diflerent temperatures (turbulence) is at the base of the boundarv laver in- 
version. 

meer has a land sea breeze even during the "trade wind" season (see Chapter 2. 
section C). [his daily shiit of circulation generating a warm, convective air mass over 
Ium the davtime followed by a cool, relatively dry environment at night renders the 
mediocre r, measurements. Katabatic flow off of the 3.1 km high Mt. Haleakala is an- 
other turbulence producing mechanism observed at the Kiher site. 

Profiling the atmosphere over Kihei during summer conditions, one expects а 
land sea breeze circulation. a boundary laver inversion, trade winds, a second moist laver 
aloft and a moderate jet stream. The shear surrounding any jet 1s a source of turbulence 
exhibited in the r, magnitudes. Shear provides the mechanical mixing necessarv for op- 
иса! turbulence, however, without the compliment of a strong temperature gradient, the 
effect on the optical regime 1s Weakened. This is best exemplified by reviewing the 
summer trade season optical measurements collected when an upper-level jet 1s travers- 
ing the site (weak temperature gradient). The 0, measurements, which are sensitive to 
turbulence in the 300 to 200 mb jet stream level, are rarely below 3 u rad. Note: 3 u rad 


I$ considered poor optical conditions. The exception to the rule for Maui occurs when 


a svnoptic-scale event crosses overhead, tightening the temperature gradient associated 
with the Jet (see Appendix D). A stronger example 15 exanuning the influence of a polar 
jet. The tighter temperature gradient and subsequently stronger winds associated with 
the polar jet are able to drive 0 measurements down to | to 2 yw rads, a very turbulent 
environment optically. An example of this was in the 10 December 1985 Haleakala data, 
where 0, was measured as 1.6 д гай. 

One of the most optically active Hawaiian cases occurs in April 1988, a winter case. 
The Kihei site experiences the passage of a low or "eddy. A detailed discussion of the 


optical and meteorological conditions throughout that experiment 1s found in Appendix 
D. 


B. OPTICAL VARIABLES CALCULATED FROM RAWINSONDE DATA 

During four of six RME experiments at Kihe1, Maur on-site rawmsonde launches 
were released. Two sessions emploved the use of a Vaisala DigiCORA MW И 
Rawinsonde Svstem. The eight April 1988 launches and ten flights recorded in Julv 1988 
began with a two-second sampling rate up to 2 km, followed bv a ten-second sampling 
rate for the remainder of the flight. All 18 vertical profile sets were utilized for цена 
four distinct C2 calculation methods. The final rating of each method was done by using 
the C} profile to calculate the two optical parameters. ry and 0,. These calculated values 
were compared to the optical measurements acquired simultaneously with the 
rawinsonde launch. 

The first two methods, the Simple Potential Temperature Differences (SC RRE 
the Bulk Gradient (BG). produced similar results. Both over estimated the turbulence 
leading to an under estimation of ry and 0, bv a factor of 10. 

Using the entire profile, the 111 and 13531 running mean techniques produced 7; and 
0, values of the same order magnitude as the sensor values. The 1 data even showed a 
general ability to match the relatively large and small magnitude variations observed in 
the sensor data. 0, calculated from the running mean methods were generaliv msensitive 
to the magnitude changes measured by the isoplanometer. However, the сајсшани 
values were within the same order of magnitude as the sensor. Overall. Ше 1355 тие 
ning mean provided turbulence parameters at least 16% closer to the sensor values than 
the 111 running mean. 

Using the Richardson number to isolate the turbulent laver from the non-turbulent 
brought mixed results. Conditioning the CZ profile by calculating turbulence for layers 


with Ri « 0.25 generated r, two to three times larger than the sensor. The DON 








provided values approximately 7% closer to the sensor. The running mean approaches 
were within 8-71 mm of the sensor. 0, calculations for SPTD, BG and running mean 
generated results such that about 40% over estimated the sensor's 0, ; 50*?o under esu- 
med. and, l0*o straddled the sensor value. 

The conditionally unstable results (Ri<2) for the SPTD Case under estimated the 
ry sensor bv 17-80?5. 

Poor correlation between the calculated and measured optical parameters revealed 
the existence of serious problems. Limitations in the current sensor system prevented 


successful C2 calculations. The limitations were broken down into three areas: 
Ei cessie Temperature Probe Response Time, 
2. Limited Vertical Resolution (aliasing) and 


Limited Temperature Resolution (Roundof?f Error). 


U- 


The limited temporal resolution of the Di;?CORA MW 11 rawinsonde profiles was 
confirmed by calculating the temperature probe time constant correction. Гог Ше 
2-second data. only 1S% of the actual Ci was measured. While the 10-second filter 
ВЕНЕЦ 5420 of the actual C$, this was deceptive. The distance between 10-second 
samples produced aliasing effects when the mean potential temperature gradient was 
large. This inability to resolve high frequencies artificially enlarged the already over es- 
timated turbulence present. In addition. this large 40-42 m separation approached the 
Kolmogorov outerscale. 

The 2-second data was not as susceptible to aliasing. Instead, the 0.1] K temperature 
resolution (rouncoff error) introduced quantization noise. Using two distinct proce- 
Mees, the standard error formula and a vertical C quantization noise simulation model. 
И оггесполе ој 6.2х 1075 А 21723 (2-зесопа Чата) апа 1.4х10-А?т-23 (10-seeond data) 
were calculated. Application of this correction confirmed that the C2 values in the initial 
2 km of data were mostly quantization noise. Because most optical turbulence is ex- 
pected near the surface, this realization was not encouraging. 

Another major problem with calculating C? from rawinsonde data was insufficient 
vertical resolution that produced a fitting problem with the running mean case. Just as 
the inability to resolve higher frequencies (aliasing) produced an over estimation of C2 
for the SPTD Case, this effect became even more apparent when the running mean 
methods encountered a 3-5 ° C inversion with a 10-second sampling rate. This enhanced 


the existing over estimation of turbulence. In contrast. turbulence calculated across the 


oF 


sharp temperature shift (inversion) captured within the 2-second « 
1988) did not appear to suffer from this difficulty. - 
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ҮП. RECOMMENDATIONS FOR FUTURE RESEARCH 


A. LINKING OPTICAL DATA AND METEOROLOGICAL EVENTS 

Quantizing the optical turbulence based on weather conditions is a non-trivial 
project. As simultaneous optical and meteorological observations accumulate, the pos- 
sibillity of plotting "weather" charts of iso — r, and iso — 0, lines nears reality. Such a 
task, however, requires at least two phases of investigation. The first deals with 
svnoptic-scale events. Major atmospheric phenomena aflecting optical propagation in- 
clude: the passage of cold and warm fronts, the polar and subtropical jets, high and low 
pressure systems, and so on. A dedicated chmatological study (spanning many vears) 
of optical parameters during each of these events, with measurements taken at one lo- 
cation, would provide a statistical sample from which ranges of values could be deter- 
mined. Being realistic, a more cost effective and time efficient approach for acquiring a 
(marginally) vald statistical sample is to design a set of experiments with multiple 
(minimum 3. preferable 5) r, and 0, sensors. The geometric shape formed by 4 of the 
sensor site locations would be a square. The sides of this square would be along the 
prevailing winds and measure a kilometer or more in length. The fifth sensor would be 
in the center of this box. Time based cross-sections would be derived from three of the 
sensors. Because the central site is conimon to both possible cross-sections, this sensor 
would be the most critical. 

The meteorological sensors needed for simultancous measurements would depend 
on the particular event being studied. For example, the passage of a cold front would 
require surface metcorological stations, rawinsondes and satellite resources. The effect 
of a polar jet would require rawinsondes, aircraft and satellite data. The main objective 
is tO pinpoint optical measurements (processed in rea] time) with the movement of the 
weather phenomena. Armed with at least three "clean" case studies of a specific weather 
event, r, and 0, ranges for various sectors of the specified event would be possible. 

The second area of investigation involves the meso- and microscale weather activ- 
ities. Though weather phenomena on these scales are often site specific, one major 
source of optical distortrons commion to all sites is the boundary laver. The influence 
of the boundary layer has been estimated as 30-50%. A simultaneous data collection 
at the surface and above the boundary laver inversion would better quantize the impor- 


tance of this atmospheric structure. Because Mt. Haleakala is over 3 km high, such a 


4] 


project could theoretically be done in Maui. However, Halcakala’s Crater experiences 
a local jet. Therefore, the measurements above the boundary laver inversion would have 
to be sampled sufficiently far from this accelerated flow. Another possible location 15 
along the central California coast where ғу values have often been recorded to be more 
than 200 mm. A dedicated study contrasting the good “seeing” conditions above and 
below the boundarv layer interface may provide a better resolution for this question. 
Another meteorological phenomena at the mesoscale is gravity waves. With the 
current method of data acquisition for the isoplanometer, tracing the influence of gravity 
waves in real-time is possible. This ability coupled with the mobility of the 


isoplanometer could be of value to gravity wave studies. 


B. CALCULATING C; FROM RAWINSONDE DATA 

Another approach to resolving the correlution of optical turbulence and 
meteorological conditions 1s to calculate key optical parameters (such ав ғ and @,) di- 
rectly from rawinsonde data. Currently, one of the main obstacles is the need for a low 
cost temperature probe and operating svstein capable of acquiring large high resoluuon 
data bases. An adequate single probe rawinsonde temperature sensor would require a 
0.01-0.02 K resolution, 100 msec response time and a ] m vertical resolution. These 
specifications are not unreasonable with current technology. Some method of detrend- 
me the potential temperature profile with a running mean still will be necessary until 
0.1-0.2 m vertical resolution is possible. 

An interesting experiment would be to launch a rawinsonde with both the above 
recommended single temperature probe and the more expensive two temperature probe 
sensor. The latter has been discussed in Olmstead (1988). A comparison of the two Сі 
profiles generated from these inethods would help to calibrate the less expensive single 
probe method. 

‘As personal computers become faster and high data-acquisition rate rawinsonde 
packages become routine. the theoretical concept of evaluating and predicting the op- 


tical turbulence over any site in real trme will evolve and become a reality. 
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APPENDIX A. CALCULATED/MEASURED OPTICAL PARAMETERS 


The optical variables r, (transverse coherence length) and 6, (1soplanatic angle) de- 
pend on an integral of C? along the vertical path. The refractive structure parameter C? 
is proportional to the temperature structure parameter C2 which can be computed from 
the rawinsonde data. Four methods for computing the temperature structure parameter 
С2 from the rawinsonde data include: 

ime Simple Potential Temperature Differences (SPTD) Equation, 


2. detrending the potential temperature differences with a 111 Running Mean (111 
px ora 13531] Running*Mean (13531 КА). 


. the Bulk Gradient Equation (BG) and 
4. conditional data sampling based on the Richardson Number (R1<0.25 and 
К: < 20). 
EC. method was tested bv calculating r, and 0, from the C$ profile. These results 
Were ultimately compared to the optical measurements simultaneously acquired over the 
Eme site. 


Table 6 and Table 7 correlates the flight numbers referenced in the Appendix А 


plots with their rawinsonde launch time (UTC). 
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Table 7. JULY 1988 KIHET, MAUI RAWINSONDE LAUNCHES 
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А. R,RESULTS 

Figure 9 through Figure 11 display the r, values calculated from the eight April 
1988 and ten July 1988 Maui rawinsonde datasets. An estimate of the r, simultaneously 
measured over the same Kihei, Maui site is also included. 


For the measured r, time series corresponding to these profiles, see Nelson (1988). 
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Ee) ©, RESULTS 

[igure 12 through Figure 14 display the 0, values calculated from the eight April 
1988 and ten July 1988 Maui rawinsonde datasets. An estimate of the 0, simultaneouslv 
measured over the same Kihei, Maui site 15 also included. 


Шетпе measured €, time series corresponding to these profiles, see Nelson (1988). 
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Tig. 13. | Measured and Calculated 0 , (Ri « 0.25) 
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APPENDIX B. MAUI, HAWAII - AUGUST 1987 


A. METEOROLOGICAL DATA 

Appendix B includes the August 1987 vertical profiles of atmospheric data acquired 
during the Relay Mirror Experiment at Kihei, Maui. The data acquisition system used 
was a VIZ WO-8000 RP+ rawinsonde system. Vertical resolution 1s approximately 35 
m. Launch site was at 20.45° latitude, 156.27° longitude. Table couples the 


rawinsonde flight numbers and launch times for all the August 1987 Kihei soundings. 


Table 8. AUGUST 1987 KIHET MAUI RAWINSONDE LAUNCHES 
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APPENDIX C. MAUI, HAWAII - DECEMBER 1987 


A. METEOROLOGICAL DATA 

Appendix C Includes vertical profiles from the December 1987 Kihei, Maui 
rawinsonde launches. The VIZ WO-8000 RP+ rawinsonde system used to acquire this 
data sampled values approximately every 35 m. Launch site for these profiles was at 
20.54° latitude, 156.27° longitude. Table 9 tabulates all the flight numbers and launch 


umes for the December 1987 Kiher soundings. 


Table 9. DECEMBER 1987 KIHEI, MAUI RAWINSONDE LAUNCHES 
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Fig. 27. Kihei, Maui Rawinsonde Data: 0950 UTC 5 December 1987 
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Fig. 36. Kihei, Maui Rawinsonde Data: 0345 UTC 9 December 1987 
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APPENDIX D. MAUI, HAWAII - APRIL 1988 


A. METEOROLOGICAL DATA 

Appendix D includes the April 1988 vertical profiles of atmospheric data acquired 
for the Relay Mirror Experiment at Kihei, Maui. The data acquisition system used was 
a Vaisala Digi CORA MW 11 rawinsonde system. Launch site was at 20.75° latitude, 
156.43" longitude. The initial 2 km were sampled every two seconds (8-10 m). Data 
acquisition above 2 km was once every ten seconds (40-45 m). Flight number and 


launch times for all the April 1988 soundings are provided in Table 10. 


Table 10. APRIL 1988 KIHEI. MAUI RAWINSONDE LAUNCHES 
ШЫ TIME (UTC) DATE (DAM MON Hm 
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2$ April 1953 
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29 April 1988 
30 April 1988 
30 April 1988 
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Fig. 38.  Kihei, Maui Rawinsonde Data: 0610 UTC 27 April 1988 
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Fig. 56.  Kihet, Maui Rawinsonde Data: 05-14 UTC 30 April 1988 
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Fig. 59.  Kihei, Maui Rawinsonde Data: 1434 UTC 30 April 1988 


100 


Aprite OO —yitlight 6 
Mixing Ratio (g/kg) 
4 6 10 


0 e 8 ра 14 
28 
26 
т 24 
E. 
с 22 за 
= 20 
Y 
E -18 
Ж Ц 
ES Е 
у: 14 
Q 12 
Š 
С) 10 
Q 7 
с 8 
` 6 





6 
4 
E 
0 


РТТ | 
100 200 300 400 500 600 700 вос 
Potential Temperature (К) 


Fig. 60. Kihei, Maui Rawinsonde Data: 1434 UTC 30 April 1988 
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B. APRIL 1988, EXPERIMENT REPORT 
The following is a general meteorological summary corresponding to the optical 
Gata measured during the 27 April 1988 (LIC) through 01 May. 1988 (UTC) Kiher, 
Maui, Hawai Experiment. 
1. 0600-1730 UTC 27 April 1988 
Extensive cloud cover prohibited. optical measurements until approximately 
1100 UTC. The 0610 UTC Kihei rawinsonde launch includes the following surface pa- 


ISmeter s: 


Ix Pressure) — 1005.9 mb 

Ши Гепшрегасигс) = 22.6%C 

RH (Relative Humidity) = 71% 

puuspeed. Wind Direction = calm (< 2m s), N/A 

Stability = unstable between surface and 0.2 km 
The stratus deck restricting optical measurements occurs between 2.2-2.7 km. and ex- 
tends horizontally to the NE 3-5" latitude according to the GOES-WEST Infrared (IR) 
Satellite Image taken within 30 minutes of balloon launch time. A temperature inversion 
of 5.1° C between 2.9 and 3.2 km caps the cloud laver. and 15 coupled with light winds 
Eu топі Пс ХЕ. 

Winds graduallv increase with height until 9.9 km (289.7 mb) where a maximum 
wind speed of 22.0 m s is reported. Wind direction remains from the eastern quadrants 
(0-180° from the North) up to 13.2 km (173.1 mb). Within the next 97m (174.3-171.1 
mb) moderate winds become calm and veer 201 degrees. The subsequent increase m 
EuuEDeed peaks at 13.7 m s (143.2 mb) before it gradually decreases to winds of less 
NU ms. Wind direction continues to be from the NW until 18.7 km (well after the 
initial signature of the tropopause) where it crosses mto the eastern quadrants for the 
remainder of the vertical sampling. The rawinsonde flight 15 terminated after 90 minutes 
ЈО seconds at 25.019 kin (24.9 mb). 

Comparing the 0546 UTC and the 1246 UTC IR Satellite Images, the dense 
Stratus deck moves from NE of the site to SW. The 1246 UTC and 1446 UTC IR Sat- 
ellite Images display very little change in the intermittent stratus deck, though r, meas- 
urements indicate a 1 cm increase over the same time interval. 

Shortly before sunrise (1548 UTC 27 April 1988) а second rawinsonde launch 


records the following surface parameters: 


Ке 550гс) = 1005.3 тЫ 
К п рсгатиге} = 1$.5°С 
RH (Ксіацуе Нипийну) = 6790 


Wind Speed, Wind Direction = сайт «2m и 
Stability = stable conditions 
The first of two inversions is found between the surface (0.083 km) and 0:257 En 
second slight temperature increase of 0.4° C (0.365-0.450 kin) is capped by a small hu- 
mudity increase of 13% between 0.459-1.012 km. The small temperature inversion 1s 
supported by 0.9 m’s increase in wind speed. Wind direction remains constant at 28° . 
Another increase їп moisture content 1s displayed between 1.270-1.737 km. The RH 
maximum of 839^ 1s reported at 1.737 km. Winds sampled thus far are light and from 
the NE, NNE. Fora thin laver, 1.782-2.275 km, the light winds back го погеплувана міни 
A sharp decrease of 52% RH (2.075-2.713 km) overlaps the variable winds. Embedded 
Within this laver 1s a 1.5° C mversion between 2.511-2.524 Кт. Except for asm 
moisture increase between 7.754-8.824 km, the remainder of the sounding displays an 
extremely drv atmosphere. At 9.520 km, winds reach a peak velocity of 17.7 ms. Ata 
slightly higher level, 10.558-12.396 km, wind speed decreases 11.9 m’s and veers 163 de- 
grees. The veering continues another 44 degrees over the next 0.943 km. The temper- 
ature gradient associated with the 217 degree wind shift 1s -20° Ci2.781 Киа s ВИ 
less than the dryv-adiabatic lapse rate). The tropopause begins around 17.3 km. 
2. 1 0000-1730 11 28 ри 1955 
Kihei rawinsonde surface measurements recorded at 0605 UTC include: 
P (Pressure $1005 в 
T (Temperature) =m C 
RH (Relative Humidity) = 72a 
Wind Speed. Wind Direccion = calir c TES Е 
Stability = unstable between surface and 0.2 km 
RH reaches a maximum of 91275 just prior to the first of two inversions. The initial in- 
Version is between 2.254-2.430 km (778.4-762.0 mb). This 2.5" C increase 15 coupled with 
an S degree veering of winds, and no change in wind speed. The second inversion Is 
fragmented into two parts. The total layer displavs a 1.1 &degee C increase over 0.3 km 
2.793-3.193 kmi). Winds associated with this stable laver increase by | пи/5 апо тесния 
degrees. As before, the inversion 1s preceded bv a slight increase (+ 1879) in RH. 
Between 5.627 km and 6.004 km light winds veer 132 degrees NEE 10 TTR 
temperature gradient associated with this large wind shift 1s СКТ ка оО 


less than the drn adia ic lapso y no 
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Another large veering in wind direction 1s between 8.951-10.830 km. This 203 
degree change is coupled with a 9.9 m s increase in wind speed and a temperature gra- 
memo. -13.2 1.879 km (2-775 Crkm less than the dry-adiabatic alpse rate). 

The jet stream occurs between 11.561-14.331 km and flows from the NNW,NW. 
E esusds of 32.1 m's at 315° are evidient at 12.493 km. The temperature gradient 
Euossuhe jet 1s -23. 7" C/2.870 km (-8.606" C;/km). The Lihue rawinsonde data recorded 
ЕШ 525 hours after the local Kiher launch displays the jet stream between 
12.354-14.580 kin (wind direction across the Jet backs from 292? to 279? ). Temperature 
gradient associated with Lihue's jet stream measurements is approximately -19.0° 
E o km(7.1643 C km). Peak wind is reported as 42.8 m’s, 286° at 13.487 km; 42.8 
185255" at 15.940 km (Note: 7j 15 5 cm during both sampling intervals 0605 and 1115 
EIC.) 

ГӘ PC 25 Apri 1988S. surface measurements are: 

mressure) = T mb 

Ши Шешрегатите) - 17.9" С 

RII (Relative n = 90% 

Speed. Wind Direction = calm (< 2m's). NA 

Stability = stable 
А 25 С inversion, as well as a 14%0 drop m RII, marks the first 84 m of the vertical 
profiles. Winds remain light and variable for the initial 2 km. A second nunor inversion 
m between 2.368 and 2.507 km displaying a temperature gradient of 0.5° C.0.438 km 
EUM C km). lhe easterly winds associated with this event increase 2.6 m's and veer 
6 degrees 

A general wind direction shift from eastern quadrants to western quadrants 
takes place around 500 mb (6.157 km), Note: 500 mb 1s often labelled a “level of non- 
divergence’. The Level of Non-Divergence is defined as a “mid-tropospheric surface 
separating major regions of horizontal convergence and divergent associated with the 
E Lerucal| structure.of the migratory cyclonic scale weather systems” (Fuschke. 
1970). In the upper half of the profile (above 300mb), RII increases at a rate of 36%o 
in 2.664 km. This secondary peak of 4975 RII is measured at 12.005 km. Nestled above 
the humidity influx is the jet stream, 12.765-14.448 km. Peak winds of 47.0 m.s are re- 
EE 15.319-13.301 km. The temperature gradient across the Jet stream 15 -12.5° 
ns km (-7.308* C' km). 


3. 0600-1730 UTC 29 April 1988 
Optical data gathering during this session was severly restricted by cloud cover 
from a synoptical scale Weather system. Unlike previous nights, the 0605 UTC sounding 
reveals the near surface winds to be from the SW. (This implies the presence of an 


anabatic flow between Krhei and Haleakala.) The rest of the surface parameters are: 


P (Pressure) = 1006.5 mb 

T (Temperature) = 24.0° C 

RH (Relative Humidity) = 75% 

Wind Speed, Wind Direction = calm (< 2m’s), N'A 

Stability = unstable between surface and 0.2 km 

The vertical humidity profile displys two peaks. The first, a gradual 13%/’km 
increase climaxing at $4%o between 1.261-1.330 km. A minor 0.1° C 0.008 m inversion 
caps this first humiditv maximum. The second peak follows a sharp 739. km gradient 
that sends the humidity soaring to 75% at 8.976 km. Most significant to this sampling, 
however, is the wind profile. As stated before. near surface winds (0.219-1.682 km) are 
hight and SW. A backing into the eastern quadrants follows and persists between 
1.739-3.104 Ка. A 0.5° C inversion interrupts the atmospheric structure Бе ЕЕ 
2.017-2.21] km. Above 5.] km, winds veer into the western quadrants where 
main for the rest of the vertical sounding. The jet stream begins at 13.046 km and ex- 
tends to 14.902 km. The maximum velocity sampled is 39.7 m’s (15721577070 
The temperature gradient accompanying the jet is -11.4° C'1.856 km (-6.142° Gen 
Using the 0546 UTC IR Satelhte Image, the site appears in the area of a comma 

cloud where one would expect overrunning and warm advection: the warm front would 
be to the south. At 1146 UTC, the site appears to straddle ше cold front. леш и 
sounding launched at 1115 UTC, records a surface pressure of approximately они 
greater than Kihei. after interpolating to 83 m (height of the Kihei site). зпасе ГИ 
approximately 2° north and 3° west of Kihei, this post-cold front characteristic 1s ex- 
pected. Associated with the approximate time of cold front crossing (and all the ex- 
pected tightening of gradients), the r, sensor displays a 20+ mm drop (!) within the same 
time interval. At 1300 UTC, r, measurements indicate a 15 mm plunge in magnitude. 
The corresponding satellite image implies the passing of additional cloud lines associated 
with smaller fractions of the larger cold front during that time period. The 1446 LTC 
image displays the passing of the most extensive cloud cover associated with the cold 
front in the comma cloud. By 1546 UTC (12 minutes prior to Kihei sonde launch 6). the 


site is boardering the area that is normally (in more northerlv latitudes) associated with 
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the dry tongue, “a protrusion of relauvelv dry air into a region of higher moisture соп- 
tent” (Huschke, 1970). In this sector, one expects cold advection and sinking arr. 
The 1558 UTC Kihei rawinsonde launch records the following surface parame- 

térs: 

P (Pressure) = 1006.6 mb 

| (Тетрегашге) = 130 С 

RH (Relative Humidity) = 86% 

Band Speed, Wind Direction = calm (< 21165), М/А 

Stability = stable 
A small inversion of 2.8? C over 57 m initials the profile. Winds remain light and vari- 
able throughout the first 5 km. Humidity measurements indicate a small peak of 72% 
at 0.561 km followed by two minor mverstons. A second peak of 71% at 1.715 km 
precedes a small 0.5° С increase between 1.768-1.815 km. A third humidity maximum 
of 62% 1s recorded at 3.159 km with vet another maximum of 69% at 7.658 km. The 
highest humidity measurement (apart from the surface) occurs at 9.307 km (73%). with 
a second strong influx of humidity at 10.151] km (72%). This moist laver immediately 
EEESedes a westerl jet stream (10.964-13.037 km). Unlike the former soundings, this 
Jets maximum wind speed is 38.6 m s... The temperature gradient across the jet stream 
ПОЛЕТО © 2.138 km (-6.36]1?^ C'km). Winds remain out of the western quadrants well 
past the initial signature of the tropopause (17 km). Not until 20 km, do the winds 
suddenlv become light and variable. 

4. 0500-1730 UTC 30 April 1988 
The 0544 UTC Kihei launch surface parameters include: 


Bee ressurc) = 1007.9 mb 

Ши perature) = 23.3° C 

ЕЕК сапе Питтак) = 79% 

Wind Speed. Wind Direction = light, W 

Stabilitv — stable 
Surface stability 1s defined bv a 0.5? C inversion in the first 6 m of the profile. An in- 
termittent laver of stratus begins at 2.082 km (96% RH). A second stable layer marked 
ру а 1.3” С increase in temperature over 0.379 km (3.430° C;km) caps the moist laver. 
05 7 km the winds pick up and veer from the NW to SE. The second of three peaks 
in the humidity is recorded at 7.113 km (S4% RII). The third humidity maximum at 
5.952 km (7520) immediately follows a laver of strong winds (33.8 m/s) from the SW. 
A weak 0.6" C 0.101 km (5.941° C, km) inversion ties this laver’s structure together. A 


Slow westerly jet stream (max wind speed of 36.8k m s at 12.972-13.021 km) occurs be- 
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tween 12.716 km and 13.228 km. This is supported by a -2.2° С/0.512% Ет (АЕ 
C km) temperature gradient. Despite a decrease 1n wind speed with height. winds remain 
W,SW up through 18 km. Around the height of the tropopause (20.5 km) winds back 
into the eastern quadrants where they remain until the sampling is terminated (23.975 
km). 

Consulting the 1115 UTC Lihue sounding, as well as the 1146 UTC satellite 
image, the overall structure of the comma cloud is now fragmented. Low level clouds 
are intermittently dispersed around the Hawaiian Islands. A moist layer between 1-2.1 
km is implied from Lihue’s sounding. Surface winds are moderate to light and from the 
east until about 6 km where thev veer and increase in magnitude. The maximum velocity 
recorded, however, is 3.0 m s from the west (12.536 km). This is hardly the исе 
of the svstem just monitored less than 45 hours ago. The r, measurements reflect this 
lack of gross turbulence with an increase from 45 mm to nearlv 70 mm sampled over the 
same 5-6 hour time period. 

The last Kiher rawinsonde launch of the experiment (1434 UTC) records the 


following surface measurements: 

P (Pressure) = 1007.4 mb 

С Гепарвева ноја 0 

RH (Relative Hunidity) 7 100?6 

Wind Speed. Wind Direction = heht. NE 

Stabilitv = stable 
The vertical profile begins with a 3.4° C inversion in the first 36 m. Surface humidity 
оГ 100% (а measurement confirmed by condensation on the telescope lenses) decreases 
20“о Іл the initial 0.312km. The light NE surface winds back to пері МЕ 
0.265km of the surface. Humidity increases to a peak of $52o at 0.6355 Km. dips about 
5*5», then increases to $7*o by 1.241 km. A 0.5° C increase occurs between [32451380 
km (4.310° C km). followed bv another inversion of 0.6° C between 2.315-2.464 km. 
Winds displav an increasing trend unul approximately 4.6 km (SW). At this same level. 
humidity plunges to values of less than 10%. The westerly jet stream 1s nestled between 
12.060-13.196 km carrving a peak velocity of 37.9 m/s betwen 12.444-12 527 kine 
tropopause occurs around 19.] km. 

5. 0600-1730 UTC 01 May 1988 
Early evening r, data indicates intermittent cloud cover. The 0546 UTC satellite 

image displavs cloud cover as a local event. The closest cyclone is NW of the site by 


approximately 10° Jatitude (1100 km). The storm system that traveled through the area 
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earlier in the week has dissolved into the high pressure ridge NE of Hawaii. The two 
available soundings (Flilo and Lihue) launched at 1100 UTC and 1115 UTC, respec- 
ШЕК correspond to relatively high +, measurements (SO mm). Both rawinsonde 
soundings indicate an absence of surface inversion. Dewpoint depression minimunis (по 
RH is recorded in these two rawinsonde datasets) occur between 0.4 and 1.6 km implying 
a moist laver within this 1.2 km region. Super-adiabatic lavpse rates are noted between 
surface-1000 mb (12.327? C; km) at Lihue. and 709-704 mb (46.867? C/km) at Hilo. 
Because there is no consistent level for these extremelv large lapse rates, thev are most 
lkelv locallv induced. Winds below 4 km are generally from the eastern quadrants 
ИС састегу winds extend to 5 km). Above this level, wind speed strengthens and 
wind direction is from the west. Jet stream velocities are recorded as W, NW between 
11.792-13.922 km over Hilo. and W.SW between 10.645-14.672 km over Lihue. Maxi- 
ши speeds of 42.8 m's at 11.755 km (Hilo) and 39.2 m's at 12.016 km (Lihue) 
ExNecorded. Doth sites measure a lapse rate of -6.7^ C/km across the jet. The two 
soundings also concur that the tropopasue caps the jet’s subsiding winds at approxi- 
ЕВ, 17.5 Кі. 

саза те тесота a sradual 30 mm decrease between 1200-1500 LTC. 
The onlv available resource over that time interval is satellite images. Apart from an 
apparent increase in local cloud cover. there is no glaring evidence of turbulent sources. 
A review of 1988 Mav 2 reveals the development of another comma cloud structure west 
атап Perhaps the sudden drop in 7, (turbulence increase) is a results of the over- 
Ene warm advection) often associated with the NE sector of this synoptic scale 


ҮСТІ. 
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APPENDIX E. MAUI, HAWAII - JULY 1988 


A. METEOROLOGICAL DATA 

Appendix E includes profiles from the July 1988 Kihei, Maui rawinsonde dataset. 
A Vaisala DigiCORA MW = 11 rawinsonde system was used as the data acquisition svs- 
tem. Launch site for this session was located at 20.75" latitude, 156.35? longitude. The 
initial 2 km were sampled every two seconds (8-10 m). Above 2 km, the sampling rate 
was once every ten seconds (40-45 m). Flight number and launch times are provided in 
Table I1. 


Table 11. JULY 1988 KIHEI. NAUI RAWINSONDE LAUNCHES 


FLECHT 


ХО ЕС) РАТЕ (ОАЖ МООИ е 


0705 7 Julx 1955 
7 Julv 1988 
o Julv 1988 


№ Š July 1988 


9 me 1955 
10 July 1955 


|| ТЕ |. 
0920 1] Лиїх 1955 
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Fig. 62. Kihei, Maui Rawinsonde Data: 0703 UIC 7 July 1988 
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Не. 63. Kihei, Maui Rawinsonde Data: 0703 UTC 7 July 1988 
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Fig. 66. Kihei, Maui Rawinsonde Data: 1537 UTC 7 July 1988 
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Fig. 67. Kihei, Maui Rawinsonde Data: 1537 UTC 7 July 1988 
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Fig. 68. Kihei, Maui Rawinsonde Data: 0536 UTC 8 July 1988 
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Fig. 69. — Kihei. Maui Rawinsonde Data: 0536 UTC 8 July 1988 
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Fig. 71. Kihei, Maui Rawinsonde Data: 1612 UTC 8 July 1988 
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Fig. 73.  Kihei, Maui Rawinsonde Data: 1612 UTC 8 July 1988 
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Fig. 74. Kihei, Maui Rawinsonde Data: 0605 UTC 9 July 1988 
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Fig. 75. Kihei, Maui Rawinsonde Data: 0605 UTC 9 July 1988 
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Fig. 77. Kihei, Maui Rawinsonde Data: 1827 UTC 9 July 1988 
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Fig. 80. Kihei, Maui Rawinsonde Data: 0556 UTC 10 July 1988 
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Fig. 81.  Kihei, Maui Rawinsonde Data: 0556 UTC 10 July 1985 
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Fig. S5. Kihei. Maui Rawinsonde Data: 1532 UTC 10 July 1988 
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Fig. 87. Kihei, Maui Rawinsonde Data: 0550 UTC 11 July 1988 
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Fig. 88.  Kihei, Maui Rawinsonde Data: 0550 UTC 11 July 1988 
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Fig. 89.  Kiliei, Maui Rawinsonde Data: 0926 UTC 11 July 1988 
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Fig. 91.  kihei, Maui Rawinsonde Data: 0926 UTC 11 July 1988 
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